16.451 Lecture 22: Deuteron Quadrupole Moment 24 Nov. 2005 


Assume that the charge distribution is an ellipsoid of revolution with symmetry axis 
along the total angular momentum vector J: 



As for the magnetic dipole moment, 
we specify the intrinsic electric 
quadrupole moment as the expectation 
value when J is maximally aligned with z: 


Qint 

= 

( e 2 ) 

rrij=J 



“Quantum geometry": cos/? 


mj_ = J 

\j\ Ji(j+ 1 ) 


If an electric field gradient is applied along the z-axis as shown, the observable 
energy will shift by an amount corresponding to the intrinsic quadrupole moment 
transformed to a coordinate system rotated through angle p to align with the z-axis 


Result: 2 


Let Q )ab be the electric quadrupole moment we measure for the ellipsoidal charge 
distribution with mj = J: 



How do we apply this to anything? 

1. A spherically symmetric state (L = 0) has Q int = 0 (e.g. deuteron S-state) 

2. Even a distorted state with J = -§• will not have an observable quadrupole moment 

3. J = 1, L = 2 is the smallest value of total angular momentum for which we can observe 
a nonzero quadrupole moment 

-> these are the quantum numbers for the deuteron D- state! 
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Deuteron intrinsic quadruple moment: 
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recall our model for the deuteron wave function: 





result for the quadrupole moment: 


Qim = — \ab\ ( r 2 


note cancellation here 


/ 


10 \ / SD 20 

= + 0.00286 ± 0.00003 bn 


— b 2 (r 2 ) 

/ DD 


A good model of the N-N interaction can fit both the magnetic moment and the 
quadrupole moment of the deuteron with the same values of a and b! -> 


Basic features of the N-N potential, via the deuteron, etc.: 
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1. Independent of the value of L, the state with intrinsic spins coupled to S = 1 
has lower energy 


-> this implies a term proportional to: 




-1/4, S= 1 


+ 3/4, S = 0 


S=1 has lower energy 


2. The deuteron quadrupole moment implies a non-central component, i.e. the potential 
is not spherically symmetric. Since the symmetry axis for Q is along J, Q > 0 
means that the matter distribution is stretched out along the J - axis: 



-> This implies a "tensor" force, proportional to 

(V r) (+ -r) 


- S 12 = -(3 


— • S 2 


Q > 0, observed 


compare: magnetic dipole-dipole 
interaction, see Griff iths problem 6.20 
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N-N interaction continued.... 
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3. There is also a spin-orbit term, as deduced from N-N scattering experiments with 
a polarized beam: 

V s _o ~ (l.s) 


(This plays a very important role also in determining the correct order of energy 
levels in nuclear spectra - more later!) 


to the N-N interaction are based on a microscopic meson 


4. Finally, all contributions 
exchange mechanism: 



Where M is a n, p, oo ... meson, etc. 

and each term has a spatial dependence 
of the form: 

-m r 

e 

V(r) = g x(spin function ) 


State of the art N-N interaction model: (not to scare anybody...) 6 


PHYSICS REPORTS (Review Section of Physics Letters) 149. No. 1 (1987) 1-89. Nonti-HoJIand, Amsterdam 

(89 page exposition of one of only ~3 state-of-the-art models of the N-N 
interaction worldwide - constantly refined and updated since first release.) 

THE BONN MESON-EXCHANGE MODEL FOR THE 
NUCLEON-NUCLEON INTERACTION* 

R. MACHLEIDT + 

Los Alamos National Laboratory, MS H850, Los Alamos , NM 87545, U.S.A. 
and 

Department of Physics, University of California. Los Angeles. CA 90024, U.S.A. 

and 

K. HOUNDE*" and Ch. ELSTER 

Insdtut fur Theoretische Kentphysik for University Bonn, NussaSee 14-16, 9-5300 Bonn , Wert Germany 
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Brace yourself - it takes 2 pages just to write all the terms down!!! 


pseudoscalar mesons: 


v J m *’ r ) = ^ m P S [(^) a > ;^+ Z K/) s <f(i 

(they use V for spin) 


(tensor operator) 


<F-6) 


scalar mesons: 


K( m «> ')— £ m.{[l - \ (5)1 n-v) + i ^ y K') + n«/)n + • s] 


(F.7) 


vector mesons: 


(spin-orbit interaction) 

K(m v . 4 = m,([l + ^ (5) ]r(m,r) - IF 2 V(m/) + y(m,r)F 2 l 

( 5 ) “ |Z 1 {m,r)/.-S - iZ(m v r)S,jj 

+ i '"*{( m * / " 1 ) 2 y (' n . r ) + IW»i) ! y (m/K • <r, 

- 4 Z,(m/)i • S - JZ(m y r)S 12 } 

f 2 

+ y ("i/K • "7 ~ aZ(mj)S i2 ) , 


more of this... 
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with 


Y(jr) - , 


Z(x) - 1 + 3 /jr + 3 /**) Y(*) , 



“Yukawa f unctions" 
and derivatives... 


l,(x)-(^)\vx + l/x } )Y(x), S u 


(/»-, • r)(tr : • r) 


(F.8) 


and 



Wc use units such that A = c = 1 (he - 197.3286 MeV fm), The use of the form factor, eq, (3,3), at each 
vertex (with n, = 1) leads to the following extended expressions: 


y.W-VK*..')- 


Az ■ 


<r-<: 


V.(-*..„r) + 


Ki~ m . 




(F.9) 


where .4,. , = A, + e, A_ .. - A„ - e. e/.-L ^ L g = 10 McV is an appropriate ch oice . 


The full NN potential is the sum of the contributions from six mesons: 


Vf.r)= 1 V a (r) 
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Parameters and Predictions: 9 


R. Machkidt et ai, The Bonn meson-exchunge model for the nucleon-nucleon interaction 


Table 14 

Meson and low-energy parameters (LEF) for the confi guration space one-boson-exchang e 


potential (OBEPR) 

deuteron properties: 


m =( MeV > 4(°eV) 

LEF Theory 


ir 14.9 138.03 1.3 

p 0.95; [6.1] 769 1.3 

t, 3 548,8 1.5 

u, 20; [0.0] 782.6 1.5 

S 2.6713 983 2.0 

a 7.7823 ‘ 550 a 2.0 

£ d (MeV). 2.2246 

P* (%) 4.81 

Q , {fin 1 ) 0.274 

(4*n) 0- 8524 


A s (fm Li ) U.K-BfrD 

D/S 0.0260 

r d (fm) 1,9691 

a, (fm) -23.751 

r t (fm) 2.662 

a, (fm) 5.423 

r, (fm) 1.759 

Impressively good agreement for ~ 10,000 experimental 
data points in assorted n-p and p-p scattering 
experiments plus deuteron observables: % 2 /d.f <~ 1.07 ! 

low energy scattering 
parameters, etc: 
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On to Finite Nuclei! 


Review: 

1. Nuclear isotope chart: (lecture 1) 

• 304 isotopes with ti > 10 9 yrs 

(age of the earth) 

• 177 have even-Z, even-N and J K = 0 + 

• 121 are even-odd and only 6 are odd-odd 

• N « Z for light nuclei and N > Z for 

heavy nuclei 

2. Elastic scattering of electrons: (lecture 7) 

• nuclei are approximately spherical 

• RMS charge radius R = 1.2 A 1/3 fm fitted to electron scattering data 

• mass ~ A , and radius ~ A 1/3 so the density M/V ~ constant for nuclei 

(«2x 10 17 kg/m 3 ), implying that nuclear matter is like an incompressible fluid 




r (fm) 

Fig. 4.3 The electric charge density of three nuclei as fitted by 
/? ch (r) = Pch/tJ + exp((r - The parameters are taken from the compila- 

tion in Barrett, R, C. & Jackson, D. F. (1977), Nuclear Sizes and Structure, 
Oxford: Clarendon Press. 




continued... 
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3. Inelastic electron scattering: (lecture 9) 

• Excited states can be identified, on a scale of a few MeV above the ground state, e.g. 


E (MeV) 


9.641 

1 

3“ 

7.6542 


9* * 

4.4389 









\Z n J T '0 + : 

L T =0 



ip/p 


4. Quantum numbers for nuclear states: 


• total angular momentum J, parity n 


• isospin, I: (lecture 13) 

for a nucleus, m x = -§• (Z-N) and I = |m x |, ie lowest energy has smallest I 


• Example: "isobaric triplet" 14 C, 14 N, 14 0: 
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Nuclear states in general: 
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• states have integer or half-integer J depending on whether A is even or odd 

• different systematics and energy level spacings for different nuclei 

• some nuclei exhibit "single particle" and others "collective" excitations 

-> different models to describe this complementary behavior 


What is the potential energy function V(r) that nuclei are eigenstates of? 


This is not an easy question! The N-N interaction is too complicated to 
solve in a many-body system: state-of-the-art can go up to A = 3! 

First approximation: a square well potential, width approx, equal to nuclear radius R: 


Assume somehow that we can treat 
the binding of neutrons and protons 
like electrons in atoms - individual 
nucleons have wave functions that are 
eigenstates of some average nuclear 
potential V(r). 

Each nucleon has a binding energy B as 
shown (E = -B) 

Kinetic energy K = (V 0 - B) 

R « 1.2 A 1/3 f m; most of the wave 
function is contained inside the well, 
so this should be approximately the 
right nuclear size... 


V(r) 


0 - 
-B 


-Vn 


<--- bound state 
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Kinetic energy of bound nucleons: 17 

Key point: once we specify the width of the well, the nucleons are conf ined, and so their 
kinetic energy is essentially determined by the uncertainty principle: 

Simple estimate: Conf ining box of side 2 f m. Ap x A x ~ h 



Px = 0 

A P* * = ^((Px-Pj 2 ) = 

-J(p 2 ) = ft/Ax 

=> (p 2 ^ = 3 (h/ Ax) 2 - 

■ 3x10' MeV 2 


K_ 

M 



0.015 




Conclusion: the motion is non relativistic; K « 15 MeV 


Now we can specify the potential parameters: 18 


What next? 

• We have a complicated system of A 
nucleons. 


• About half of them are protons, so a 
repulsive (+ve energy) term has to 
be added to the square well to 
account for this (~ few MeV) 

How to connect this model to something 
observable? 


Independent particle model: 

• Assume independent particle motion in some average nuclear potential V(r) as shown. 

• Then we can fill the eigenstates of the potential to maximum occupancy to form a 
nucleus, as is done with electrons in atoms (to 1 st order...) 
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Connection to average nuclear properties: 19 

• The binding energy of each nucleon, in our model, is a few MeV. 

• The potential energy of a bound nucleon is negative, by - 0.3% of its rest mass energy, 

which therefore has to show up as a decrease in its mass. 

• For A nucleons, the total binding energy is: 



mass of nucleus, M 


The average binding energy per nucleon, B/A, can be determined from mass data 
and used to refine a model for V(r); it ranges systematically from about 1-9 MeV 
as a function of mass number for the stable isotopes. 


Reference: F<£H ch. 16 


Atomic Mass Units: 
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• By convention, we set the mass of the carbon-12 atom as a standard. 

• Denote atomic masses with a "script" M, measured in atomic mass units, U 

M ( 12 C) = 12.0000000000 U (exact!) 1 U = 931.494 MeV (expt) 


Calculation for carbon-12: 


m p = 938.2 MeV ' 
m n = 939.6 MeV y 
m e = 0.511 MeV , 


6x7(11; = 11,269.8 MeV 

Z 

12 U = 11,178.0 MeV 


1 


B ( 1 |C) = = 91-8 MeV 


Binding energy per nucleon in 12 C: B/A = 7.8 MeV; 
Contrast to the deuteron 2 H: B/A = 1.1 MeV 
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The famous Binding Energy per Nucleon curve: 


21 



Implications of the B/A curve: 


Greater binding energy implies lower mass, greater stability. 

Energy is released when configurations of nucleons change to populate the 

larger B/A region -> nuclear energy generation, e.g. 

Fusion reactions at small A release substantial energy because the B/A curve rises 
faster than a straight line at small A: 


Binding energy 
of products is greater 
* than the sum of the 
binding energies 
of the initial species. 


Fission reactions at large A release energy because the products have greater 
binding energy per nucleon than the initial species: 

235 U + n => 100 X + 133 Y + 3n + 200 MeV 


distribution of final products 


2 H+ 2 H => 3 He + n + 3.27 MeV 
=> 3 H + p + 4.03 MeV 
2 H + 3 H => 4 He + n + 17.6 MeV 
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A semi-empirical model for nuclear binding energies: 


1. Volume and Surface terms: 

First consider a 1-dimensional row of nucleons with interaction energy per pair = e 


— 0 -#- 0 -# 0 ( 
V. 


‘ 0-*0 


"V 

total: A 


A 

“1 

Is - A = 2s A - L 

\ 1 

f 

1 

1 1 

1 


each has 2 neighbors 


correction 
for the ends 



Approximately constant, with end 
effects relatively smaller at large A. 


By analogy, for a 3-d nucleus, there should be both volume and surface terms 
with the opposite sign, the surface nucleons having less binding energy: 

B = a v A - a s A 2/3 => ^ = a v - a s A~ 1/3 


2 . Coulomb term: 
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for a uniform sphere, 

r q(r) dq 3 (Zef 
Coul J 4ns 0 r 5 4 ns 0 R 

total charge + Ze 


This effect increases the total energy and so decreases the binding energy. 
Simple model: A B = - a c Z 2 A~ 1/3 


But this is not quite right, because in a sense it includes the Coulomb self energy 
of a single proton by accounting for the integral from 0 to r p ~ 0.8 f m. The nucleus 
has fuzzy edges anyway, so we will have to fit the coefficient a c to mass data. 


Solution: let AB scale as the number of proton pairs and include a term: 

AB = -a c Z(Z -1) A“ 1/3 => ^ = - a c Z (Z -1) A“ 4/3 
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3. Symmetry Term: 
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So far, our formula doesn't account for the tendency for light nuclei to have Z = N. 
The nuclear binding energy ultimately results from f illing allowed energy levels in a 
potential well V(r). The most efficient way to fill these levels is with Z = N: 


Simplest model: identical nucleons as a Fermi gas, i.e. noninteracting spin- -§■ 
particles in a box. Two can occupy each energy level. The level spacing ~ 1/A. 

A mismatch between Z and N costs an energy price of AE at fixed A as shown. 


energy 


o i 1- 


dAt 

\ r 

t 1 

At 

ri 

1 — + U Z. - IN 

neutrons 

. U states full below 

protons f , 

r J here.. 





A B = -a A (Z-N) 2 A' 1 = -a A (A- 2Z) 2 A -1 


4. Pairing Term: 


Finally, recall from slide 1 that for the case of even A, there are 177 stable 
nuclei with Z and N both even, and only 6 with Z and N both odd. Why? 

-> Conf igurations for which protons and neutrons separately can form pairs 
must be much more stable. All the even-even cases have = 0 + , implying 
that neutrons and protons have lower energy when paired to total angular 
momentum zero. 

Solution: add an empirical pairing term to the binding energy formula: 


II 

III 

! 

od 

< 

”+l” 

0 

a A' 374 

pu.li 

_-i_ 

p 


with +1 for even-even, 0 for even-odd, and -1 for odd-odd 

Full expression: 

B(Z, A) = a v A-a s A 2/3 -a c Z (Z-l) A“ 1/3 — a A (A-2Z) 2 A -1 + 8 
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Fitting of coeff icients to data: 27 



Mass number A 
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